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Abstract. IGR J19140+0951 was discovered by INTEGRAL in 2003 in the 4-100 keV band. Observations with 
INTEGRAL and RXTE provide a tentative identification as a high-mass X-ray binary (HMXB) with a neutron 
star as accretor. However, an optical counterpart was thus far not established, nor was the presence of a pulsar 
which is commonly observed in HMXBs. We observed IGR J19140+0951 with Chandra and find the source to 
be active at a similar flux as previous measurements. The lightcurve shows a marginally significant oscillation at 
6.5 ks which requires confirmation. We determine a sub-arcsecond position from the Chandra data and identify 
the heavily reddened optical counterpart 2MASS 19140422+0952577 in the 2MASS catalog. Optical follow-up 
observations with the William Herschel Telescope at La Palma exhibit a continuum spectrum coming out of 
extinction above 7000 A without strong absorption or emission features. V, I and K a band photometry point to 
an optical counterpart extincted by Av = 11 ± 2. The extinction is consistent with the interstellar value. None 
of the data reject the suspicion that IGR J19140+0951 is an HMXB with additional circumstellar obscuration 
around the accretor. 
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1. Introduction 

IGR J19140+0951 was discovered with INTEGRAL in 
March 2003 (Hannikainen et al. 2003) during observations 
of the nearby (1?1) microquasar GRS 1915+105. In obser- 
vations that continued until May the source was detected 
70% of the time above a threshold of 9 to 10 mCrab (20-40 
keV); the brightest flux was measured during a flare peak- 
ing at 70 mCrab (Hannikainen et al. 2004). The position is 
coincident with that of EXO 1912+098 which was discov- 
ered in archival EXOSAT data by Lu et al. (1997). RXTE 
had made pointed observations of the source prior to the 
report of the discovery. A public target-of-opportunity ob- 
servation was triggered in April 2002 by a detection of the 
source in the BeppoSAX Wide Field Cameras (WFCs; In 
't Zand et al. 2004). RXTE observed the source again in 
2003 after the INTEGRAL detection, for 2.8 ks. Swank & 
Markwardt (2003) find the source to vary on time scales 
longer than 100 s and peak at 10 mCrab (2-10 keV). 
is fairly high at 6 x 10 22 cm -2 , while the continuum can 
be described by a power law spectrum with an Fe-K line 
of equivalent width 500 eV. 

The persistent nature of the source, as suggested from 
the EXOSAT and WFC detections, was confirmed with 
the RXTE All-Sky Monitor. Corbet et al. (2004) discov- 



ered a sinusoidal modulation in the flux with a modula- 
tion depth of about 75% and a period of 13.558+0.004 d 
present since the start of observations in 1996. The per- 
sistence supports its identification with the orbital period 
of IGR J19140+0951. This orbital period is highly sug- 
gestive of a high-mass X-ray binary (HMXB). An HMXB 
is a binary consisting of a compact object and a normal 
star heavier than a few solar masses that is transferring 
mass to the compact object through a wind or Roche lobe 
overflow resulting in large amounts of X-ray emission. 

Rodriguez et al. (2005) made a comprehensive analy- 
sis of the March-May 2003 INTEGRAL data (1.3 Msec 
of exposure at energies above 4 keV) and 2002-2004 
RXTE observations (13 ks above 3 keV). They identify 
4 states based on the observed 20-40 keV photon flux, 
and are able to model the spectrum of each state with 
a Comptonized component of plasma temperature 10 to 
20 keV which in all but the brightest state is comple- 
mented with a faint blackbody component. Alternatively, 
the Comptonized component can be replaced by a cutoff 
power law with a photon index between 1.4 and 2.4. In 
almost all observations an Fe-K line is required. The data 
do not allow an accurate measurement of the line width, 
although the INTEGRAL data for the middle two states 
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Fig. 2. (Top) Chandra spectrum (crosses) as derived from 
the PSF with the best fit power law model (histogram). 
(Bottom) residuals with respect to the best fit model. 



Fig. 1. Map showing source detections around IGR 
J19140+0951 . There are two positions for EXO 
1912+097. The dashed circle refers to the coordinates as 
provided by Simbad, the solid circle to those inferred from 
Lu et al. (1997). The INTEGRAL data (JEM-X and IBIS 
combined) are from Cabanac et al. (2004). The asterisk 
designates the Chandra position. The WFC error circle is 
from In 't Zand et al. (2004). Circle radii are for a 90% 
confidence level. 



suggest a broad line of width 0.4-1.2 keV or more (ler). 
The source is detected up to ~ 100 keV and the unab- 
sorbed flux at maximum is 2.1 x 10 -9 erg cm _2 s _1 (1-20 

(20-200 keV). The RXTE 



keV) and l.OxlO -9 erg cm _2 s _1 
data show a variable TVh between 3 and lOxlO 22 cm~ 2 , 
bringing IGR J19140+0951 marginally in the range of the 
'obscured INTEGRAL sources' (for a recent review, see 
Kuulkers 2005). The Galactic Nu measured by Dickey & 
Lockman (1990) within 1° from the source is between 1.3 
and 1.9 x 10 22 cm -2 ; the data point closest to the source 
(at 0?24) is 1.8 x 10 22 cur 2 . Rodriguez et al. find that 
the spectroscopic evidence favors a neutron star over a 
black hole for the nature of the accretor: in the 20-200 
keV luminosity versus 1-20 keV luminosity diagram (first 
presented by Barret et al. 1996), the source is in the X-ray 
burster domain as long as the distance is below 10 kpc. 

Thus, all evidence points to IGR J19140+0951 being 
an HMXB with strong and variable local absorption and 
a neutron star as accretor. Prior to the INTEGRAL dis- 
covery we had already proposed a Chandra observation of 
EXO 1912+098, based on the WFC detection, to search 
for possible long X-ray pulsations in a faint source and 
determine an accurate position to search for the optical 
counterpart. The data were obtained in 2004. They now 
provide an opportunity to make a conclusive identification 
of IGR J19140+0951 as an HMXB. Here we present the 
results of this observation and of the optical follow up. 



2. Chandra observation 

Chandra observed IGR J19140+0951 on May 11th, 2004, 
starting at 17:53:30 UT, with the ACIS-S CCD array 
(Garmire et al. 2003) in the focal plane and no grat- 
ing. The CCD frame time of the data is 3.2 s, the ex- 
posure time 20.10 ks. The binary orbital phase cover- 
age is 0.833-0.850 with an uncertainty of 0.055 based on 
the ephemeris of Corbet et al. (2004). Thus, the obser- 
vation was well ahead of the maximum orbital modula- 
tion. The source is clearly detected. It is the only source 
brighter than 0.0003 c s" 1 inside the INTEGRAL error 
circle (Hannikainen et al. 2004) . The image shows a piled- 
up source with readout trails but no obvious hole at the 
center of the point spread function. The average posi- 
tion of all photons within 10" is a 2 ooo.o = 19 h 14 m 4?232, 
<5 2 ooo.o = +9°52'58'.'29 (l H = 44?30, b 11 = -0?47) with 
a nominal uncertainty of 0'.'6. This position is consistent 
with that of EXO 1912+098, IGR J19140+0951 and the 
BeppoSAX/WFC position, see Fig. □ 

We reprocessed the level- 1 data with CI AO version 
3.2.1 to resurrect afterglow events, using among other 
tools acis_run_hotpix. This procedure adds 6.9% to the 
number of photons within 2'.'0 of the source. We extracted 
a source spectrum from the photons in this region with 
the CIAO tool psextract; a background spectrum was 
extracted from a circular region with a radius of 1' on the 
same CCD and centered 70" from the source centroid. The 
spectral analysis was performed with XSPEC version 11.3.1 
(Arnaud 1996). We employed the pileup model in XSPEC 
which was specifically designed for the analysis of ACIS 
data (Davis 2001), leaving free the morphing parameter a 
and the fraction of the point-spread function (PSF) that 
is piled up. The spectrum was binned so that at least 15 
photons are contained per bin, thus allowing the use of the 
X 2 statistic as a goodness of fit estimator. The resulting 
spectrum can be well fitted with an absorbed power-law 
(xf, = 1.005 for 130 degrees of freedom). The photon in- 
dex is r = 1.1 ± 0.8 (all errors quoted in this paper are for 
a 90% confidence level) and A H = (1.0 + 0.3) x 10 23 cm" 2 . 
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Fig. 3. Light curves of IGR J19140+0951 during the 
Chandra observation in full bandpass. (Top) From pho- 
tons in the readout trail (resolution 1000 s) (Bottom) From 
photons within an annulus with an inner radius 2"and 
outer radius 5" (resolution 500 s). 



These values are similar to those found by Rodriguez et 
al. (2005). The absorption column is considerably larger 
than the interstellar value. Figure |21 shows the spectrum 
and model fit. Very few photons are detected below 1 keV 
and they were not taken into account when fitting models, 
neither were photons above 10 keV. The unabsorbed 1-10 
keV flux is 3.9^J'g x 10 -11 erg cm^s' 1 , but this number is 
susceptible to systematic errors because of the uncertainty 
of the pileup correction. Therefore, we verified it with 
the spectrum of the readout trail. This spectrum was ex- 
tracted with the CIAO tool acisreadcorr. A background 
spectrum was extracted from a strip on the detector 50" 
to the north. The spectral shape is consistent with the 
above result: T = 0.7 ± 0.7, Ah = (0.8 ± 0.3) x 10 23 cm" 2 
(xt = 0.96 for 51 degrees of freedom). If we fix the 
spectral shape to the more accurately determined former 
shape and fit the flux, we find (1.5 ± 0.2) x lO" 10 (1- 
10 keV) or (4.0 ± 1.0) x 10~ 10 erg cm^s" 1 (1-20 keV). 
This is fainter than in all INTEGRAL detections ana- 
lyzed by Rodriguez et al. (2005) but 50% brighter than 
in the brightest pointed RXTE observation. We find no 
evidence for an Fe-K line at 6.4 keV. Adding a Gaussian 
line at 6.4 keV improves the fit by A\ 2 — 4.5 in the PSF 



spectrum which is not significant (10% chance probabil- 
ity). The 90%-confidence upper flux limit on a narrow 
line is 7 x 10~ 5 phot s _1 cm -2 ; on a Gaussian line with 
a = 0.4 keV it is 3.7 x 10~ 4 phot s^cixT 2 . The limits 
on the equivalent widths are 116 eV and 839 eV respec- 
tively. These values are consistent with the INTEGRAL 
and RXTE detections for a broad line, but inconsistent for 
a narrow line. Thus, a comparison of the Chandra with the 
INTEGRAL/RXTE observations rules out a narrow line, 
if the line is present with a similar equivalent width. 

We generated light curves from the trailed photons at 
various time resolutions between 40 /is and 1000 s. The 
high-resolution data do not show interesting features. We 
show in Fig. |3] (top) the time profile at 1000 s resolu- 
tion. The source is obviously variable. There is even a 
hint of an oscillation with a period of about 6.5 ks. We 
also generated a light curve from photons in an annulus 
around the PSF (Fig. bottom), thus minimizing pile up 
effects. The modulation is also visible here. We checked 
RXTE All-Sky Monitor (ASM) data for a confirmation. 
The ASM data are marginally qualified for that because 
the detection threshold is about 40% , and the time resolu- 
tion is high enough (90 s) to resolve a 6.5 ks modulation. 
We tested the full data set as well as orbital phases at 
which the flux is highest and the three energy channels 
separately. The times were corrected for the Earth's orbit 
around the Sun, but not for the satellite or binary orbit 
which are not important for a possible 6.5 ks pulsation. We 
are unable to obtain an unambiguous confirmation from 
the ASM data. This negative ASM result implies that the 

6.5 ks modulation probably is not related to a neutron star 
spin. However, additional confirmation may be worthwhile 
since modulation patterns may change (e.g., Patel et al. 
2004). 

3. Optical and near-infrared observations 

3.1. Identification of counterpart 

We cross checked the Chandra position with the USNO 
B1.0 (Monet et al. 2003) and 2MASS (Skrutskie et al. 
1997) catalogs and found a 2MASS object consistent with 
the Chandra position (2MASS 19140422+0952577). It is 
fairly close (4'.'0) to another bright object to the south. 
We acquired a near-infrared image of the field in the K s 
band under non-photometric conditions using the 3.6 m 
ESO NTT with SOFI, see Fig.H It clearly reveals the two 
objects near the Chandra position. After calibrating the 
plate scale through 9 isolated stars that were identified 
in the 2MASS catalog, the northern object was localized 
at a 20 oo.o = 19 h 14 m 4?231, 5 200 o.o = +9°52'58'.'35 which is 
0'.'06 from the Chandra position of the X-ray source and 
0'.'6 from the 2MASS position of the optical counterpart. 
The southern object is present in the USNO B1.0 catalog 
while the northern is not, testifying to the redness of the 
optical counterpart. The 2MASS magnitudes are K s = 

7.06 ± 0.20 for the northern and 6.274 ± 0.026 for the 
southern object. 
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Fig. 4. (Left) i*T s -band image with the O'.'6-radius Chandra error circle. The stars designated with white crosses were 
used for the astrometric solution through cross correlation with the 2MASS catalog. This image was taken with SOFI 
on the NTT at the ESO on May 18, 2005, under non-photometric conditions. The seeing was 1'.'3 (FWHM). (Right) 
15"xl5"zoomed in image with an adapted grey scale. 



3.2. Spectroscopy 

We obtained two spectra of 2MASS 19140422+0952577 on 
August 12, 2004, each with an exposure time of 1500 sec- 
onds, using the ISIS spectrograph mounted on the 4.2 m 
William Herschel Telescope located at the Roque de Los 
Muchachos Observatory, La Palma, Spain. We used the 
R316R grating with the central wavelength set to 7500 A 
and a slit width of 1" in combination with the MARCONI2 
CCD. This setup yields a mean dispersion of 0.83 A per 
pixel. The source was observed at an airmass of ^1.1. It, as 
well as the southern source in the 2MASS image, is not re- 
solved into multiple sources with a seeing of 1'.'2 (FWHM). 
Due to the high interstellar extinction towards the source 
we are not able to detect the source spectrum in the blue 
arm of the ISIS spectrograph. In the red arm of the spec- 
trograph we detect the source redward of w 7000 A, sec 
Fig. El We find weak evidence for absorption lines from the 
Paschen series. However, due to the relatively low signal- 
to-noise ratio, the presence of telluric lines at or close to 
the wavelength of some of the Paschen lines, and fring- 
ing in the red part of the CCD for which no calibration 
measurements are available, the continuum level is too 
uncertain to determine. We estimate a conservative upper 
limit to the equivalent width of any narrow line of 1.6 A. 
The spectroscopic data have too narrow a bandpass to 
confidently determine the spectral type. The presence of 
Paschen absorption lines, if confirmed, would point to an 
A, B or late O-type companion star. 



3.3. Photometry 

The 2MASS magnitudes for the optical counterpart are 
J = 8.55 ± 0.05, H = 7.67 ± 0.16 and K s = 7.06 ± 0.20. 
However, this is based on aperture photometry without 
taking into account the blending by the southern source. 




7500 8000 
Wavelength (A) 



Fig. 5. Optical spectrum of 2MASS 19140422+0952577, 
taken with the red arm of the ISIS spectrograph on the 
WHT at La Palma. The background was subtracted, but 
no flux calibration was performed because no observa- 
tions of standard stars were performed. The instrument 
response is fairly flat between 5500 and 8500 A after which 
it drops by 60% at 8900 A and by 98% at 9500 A. Thus, 
this spectrum illustrates the strong reddening. 



Therefore, the brightnesses are overestimated. This is 
shown through a measurement of K s using the new image 
taken on May 18, 2005, albeit under non-photometric cir- 
cumstances (Fig.0J). After calibrating the magnitude scale 
with 4 isolated stars that are also present in the 2MASS 
catalog and have the same magnitude differences among 
them, wc find K s — 8.69 ± 0.05 which is 1.6 mag fainter 
than in the 2MASS catalog. 

We obtained additional images in the V (180 s to- 
tal exposure) and / band (240 s) with the 1 m Henrietta 
Swope telescope at Las Campanas Observatory on August 
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7, 2005. The seeing was 1'.'3 at an airmass of ~1.3. The 
counterpart was detected in / at 13.0±0.1 mag, while it 
was not detected in V. A limit of V > 18.8 mag was de- 
termined from the magnitude of the faintest unambiguous 
star visible in the field. 

We also checked infrared data taken with the 
Midcourse Space Experiment (MSX; Mill et al. 1994) and 
find a 0.281 ±0.013 Jy object in the most sensitive 8.3 /xm 
band at a 20 oo.o = 19 h 14 m 4!16, 62000.0 = +9°52'55"4 with 
a 1-sigma uncertainty of 0'.'8. This is 2'/5 from the optical 
counterpart (or 3a) and 1'.'6 (2a) from the nearby south- 
ern object. Given the uncertainty we cannot rule out an 
association with the unrelated southern object. 

We sought confirmation on the suspected early spectral 
type by fitting photometric models for main sequence and 
supergiant stars (Cox 2000), using extinction laws from 
Schlegel et al. (1994), to the VIK S photometry. The re- 
sults are ambiguous, see Fig. The best fit models are 
late-type stars (G3I or V), but the fits are not accept- 
able. Photometric fits with earlier spectral types are worse. 
Possibly the non-simultaneity of the K s measurement with 
the V and / measurements plays a role. 

The extinction leaves an unambiguous imprint on the 
photometry, and is fairly independent of spectral type. 
The values for Ay range between 9 and 12 and are, based 
on the translations determined by Predehl & Schmitt 
(1995), consistent with N H = (1.9 ± 0.3) x 10 22 cm~ 2 . 
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4. Discussion 

Our observations result in the secure identification of the 
optical counterpart to one of the 'obscured INTEGRAL 
sources' (e.g., Kuulkers 2005). Like the counterpart of the 
prototypical IGR J16318-4848 (Filliatre & Chaty 2004), 
it is a peculiar highly reddened object with bright NIR 
magnitudes whose spectral type still needs to be resolved. 

The most likely scenario is that the optical counter- 
part is an early-type star in an HMXB. In analogy to 
other HMXBs, the X-ray persistence, the orbital mod- 
ulation and the lack of a short pulse period suggests a 
supergiant donor star. An alternative to an HMXB sce- 
nario would be a low-mass X-ray binary (LMXB) where a 
Roche-lobe filling sub-solar mass star is donating matter 
to a compact object via an accretion disk. This is unlikely 
on two grounds. First, the orbital period is much longer 
than most LMXBs (c.f., Liu et al. 2001). Second, optical 
counterparts to LMXBs typically have Mr s between -2 
and 2 (e.g., Wachter et al. 2005). This would imply for 
IGR J19140+0951 a distance between 30 and 90 pc and a 
maximum 1-200 keV luminosity of 3 x 10 33 erg s" 1 . This 
is inconsistent with its persistently active behavior: accre- 
tion disks become unstable at such low luminosities for 
the orbital period measured for IGR J19140+0951 (e.g., 
Van Paradijs 1996). 

Support for the HMXB scenario also comes from the 
position: IGR J19140+0951 is in the direction of the tan- 
gent to the Sagittarius arm. HMXBs are expected to be 
coincident with spiral arms because that is where young 



Fig. 6. Photometry in V, I and K s (crosses; horizontal 
bars indicate band width and vertical bars magnitude er- 
rors), fitted model (solid line) and fitted model before ex- 
tinction (dashed line). The top panel shows the best fit 
with spectral type G3I and Ay = 9.2; the bottom panel 
shows, as example, the fit with a 091 model leaving free 
only Ay. Both models are unacceptable, but xt is two 
times larger for the 091 model. The models were taken 
from Cox (2000) and the extinction laws from Schlegel et 
al. (1994). 



systems reside. If the association with the Sagittarius arm 
is true, the distance is of the order of 2 to 6 kpc. The im- 
plied 1-20 keV luminosity during the Chandra observation 
is roughly 10 35 erg s , a common value for an HMXB. 

There is no evidence for obscuration of the optical 
counterpart beyond the interstellar values provided by 
Dickey & Lockman (1990), while there is such evidence 
for the accretor through X-ray observations. The orbital 
modulation of the X-ray flux (Corbet et al. 2004) sug- 
gests either the presence of an inclined circumstellar disk 
formed by the donor or an eccentric orbit, or both. 

Further understanding of IGR J19140+0951 would 
benefit from a longer Chandra or XMM-Newton observa- 
tion, to study the variability (in particular possible pul- 
sations) without data gaps due to earth occultations, and 
from a calibrated infrared spectrum in the I to L bands, 
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in order to measure spectral lines and achieve a reliable 
spectral classification. 
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